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FIELD OF THE INVENTION 

The present invention relates to a modelling method allowing to predict, as a 
function of time, the detailed composition of fluids produced by an underground 
reservoir under production, and notably a hydrocarbon reservoir under production. 

5 Such a model allows reservoir engineers to establish detailed compositional profiles 

that are in particular necessary for predictive dimensioning and management of surface 
installations, such as separators, treating plants, transport lines, etc., and therefore useful 
in surface and process engineering. 

BACKGROUND OF THE INVENTION 

10 Modelling of flows in a petroleum reservoir is essentially based on the application, 

to the previously gridded reservoir (or to a portion thereof), of the well-known Darcy's 
law describing the flow of fluids in porous media, of the law of mass conservation in 
each volume unit, of the thermodynamic relations that govern the evolution of the phase 
parameters of the fluids such as viscosity, density, on initial conditions, on structural 

15 closure boundary conditions and on bottomhole conditions. 

The « Black Oil » model, referred to as B.O. hereafter, is one of the most commonly 
used models in petroleum simulation. It allows to describe a compressible three- 
dimensional and three-phase (water-oil-gas) flow. The petroleum effluents involved in 
this model are generally described by a water constituent and by two constituents for the 
20 reservoir fluid, the term constituent covering here the notion of component (such as H 2 0 
for water) and the notion of pseudo-component (group of components). The constituents 
involved in this model are three in number: a water constituent (E), a heavy 


hydrocarbon constituent (L) and a light hydrocarbon constituent (V). In a B.O. type 
model referred to as « strict », constituent (E) is present only in the water phase, 
constituent (L) is present only in the liquid hydrocarbon phase (referred to as oil or 
condensate), and constituent (V) is distributed among the liquid and vapour 
hydrocarbon phases (phase referred to as gas). Although the use of B.O. models is not 
recommended in certain condensate gas cases, it is however applicable to a large 
number of industrial cases. 

Another well-known simulation model, referred to as « compositional » model, is 
also used, wherein the hydrocarbon fluids are represented by a larger number of 
components, at least three, often more. Modelling the flows of these more detailed 
fluids leads to very long computing times (much longer than those required for B.O. 
type modelling) as a result of the larger number of constituents, but also because it is 
often necessary to reduce the size of the grid cells to limit numerical errors and 
consequently to increase the number of cells. 

For practical reasons, the fluids in place are described as consisting of a number of 
components or pseudo-components that is much smaller than the real number of 
components, so that the modelling computations can be carried out within a reasonable 
period of time. A composition reduced to some 5 or 10 pseudo-components is generally 
sufficient to represent the behaviour of the fluids in the reservoir. 

Patent application WO-99/42,937 and the paper by C. Leibovici and J. Barker « A 
Method for Delumping the Results of a Compositional Reservoir Simulation », SPE 
49068, presented at the SPE Annual Technical Conference and Exhibition New Orleans, 
27-30 September 1998, describe a method for predicting the evolution of the detailed 
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composition in time, from computations carried out in a compositional type simulation 
of fluids described by a certain reduced number of pseudo -components (principle of 
« lumped » representation obtained by means of a « lumping » operation), the number 
of components being at least three. The method thus allows to predict the results that 
5 would have been obtained with a reservoir simulation using a finely detailed model 
where the fluids are represented by a larger number of components. This operation is 
well-known to the man skilled in the art as « delumping ». 

The principle of the prior delumping stage consists in calculating coefficient AD 0 
and the n coefficients AD P (i.e. n+1 coefficients, n being the number of parameters of 
10 the equation of state) of a known general equation, previously published in a paper by 
C.F. Leibovici, E.H. Stenby, K. Knudsen, « A Consistent Procedure for Pseudo- 
Component Delumping », Fluid Phase Equilibria, 1996, 117, 225-232 : 

(1) Ln{k) = &D (i + f u /±D p Il p! 

where the n pi are fixed characterization parameters of constituent i in the equation of 
state for a given thermodynamic representation, from the equilibrium constants k, of 

15 each constituent of the lumped thermodynamic representation computed during 
compositional simulation in each grid cell and at each time interval. If N rg is the number 
of components of the lumped thermodynamic representation, we thus have N rg equations 
to determine n+1 coefficients. A necessary condition is therefore that N rg is at least 
equal to n+1. For Peng-Robinsons' s two-parameter equation of state, a lumped 

20 thermodynamic representation with at least three components is therefore required. 
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Once coefficient AD 0 and the n coefficients AD P calculated, they are used for 
calculating the equilibrium constants of the components of the detailed thermodynamic 
representation (N rd components) by applying Equation (1) to the N rd components with 
their own fixed characterization parameters in the detailed thermodynamic 
5 representation. 

By using a) the equilibrium constants thus determined for the detailed 
thermodynamic representation, b) the flows between each grid cell and in the wells, c) 
the vapour fraction in each grid cell from the lumped compositional simulation, and d) 
the global detailed composition in each grid cell and in the injection wells at the 
10 beginning of each time interval, the detailed composition of each hydrocarbon phase at 
the time interval t and the global detailed composition of each cell at the next time 
interval (t+1) are then estimated in each cell. 

One of the advantages of this method is that it is not necessary, in the delumping 
stage at each time interval, to solve the equation of state, whether for the lumped 
15 representation or for the detailed representation, which allows to save computing time. 
One drawback of the method is that it is not applicable to B.O. type simulations since 
convenient equations of state have at least two parameters. 

Patent application WO-98/5,71 0,726 describes a method for predicting the 
evolution of the detailed composition in time from the flow computations carried out in 
20 a B.O. type simulation where the hydrocarbon phases are described by only two 
components (L) and (V). The drawback of this method is that it requires the use of the 
equation of state of the detailed representation at each time interval, and it is therefore 
time-consuming during the delumping stage. 
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SUMMARY OF THE INVENTION 

The method according to the invention combines a predictive simulation of the 
behaviour of the reservoir by means of B.O. type modelling, which is advantageous in 
that it requires no equation of state, and a delumping stage wherein no solution of an 
5 equation of state is required at any time interval. 

In B.O. type modelling, the thermodynamic parameters of the fluids (such as 
viscosity, volume factor, density, gas-oil dissolution ratio, etc.) are entered, among other 
parameters, in an input data file of the simulator (in form of charts, and/or by 
correlation, as a function of the pressure, of the temperature if it varies) as well as, if 
10 need be (for the extended B.O. type model), an additional parameter keeping a memory 
of the composition of the gas such as, for example, the density of the gas. 

The basic idea here is to simply enter in the input data the (n+1) additional columns 
or correlations corresponding to differences AD P , for example by including them in a 
file incidental to the master file used by the simulation. It is also possible to enter 2 x 
15 (n+1) additional columns or correlations corresponding to parameters D p of each phase. 
A practical condition for data to be entered in form of charts is that the functional 
dependencies introduced are sufficiently continuous. 

Computation of the equilibrium constants by the model is thus no longer required. 

The simulation method according to the invention allows to predict, as a function of 
20 time and in at least one thermodynamic zone, a detailed composition of a fluid produced 
by a reservoir and, more particularly, a detailed composition of a fluid contained in and 
produced by an oil reservoir in which one or more production wells are installed. 
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It is characterized in that it comprises the following stages : 

- discretizing the reservoir by means of a grid, each grid cell containing one or more 
phases, including at least one non-aqueous phase ; 

- determining the variation of thermodynamic parameters of the non-aqueous phases 
5 necessary for « Black Oil » (B.O.) type modelling during stages of a thermodynamic 

path followed by the fluids in the reservoir (such as, for example, the viscosity of the 
phases, the gas-oil dissolution ratio, the volume factors of the oil or the compressibility 
of the gas, etc.) ; 

defining the fluids by a detailed representation with N ri components and/or pseudo- 
10 components ; 

- determining at least one equation of state with n parameters allowing to simulate, at 
least during the input data preparation stage (a simulation that is strictly necessary only 
during this input data preparation stage, and therefore not repetitively necessary during 
or after the B.O. type simulation), the thermodynamic behaviour of the fluids during the 

15 stages of a thermodynamic path, allowing to generate, per thermodynamic range or 
zone, n+1 additional functional relations (in form of charts or of correlations) ; 

- converting the thermodynamic behaviour of each non-aqueous phase into input data 
suited for a B.O. type thermodynamic representation, this input data being completed by 
the additional input data suited for the delumping operation ; 

20 - carrying out B.O. type modelling, which allows to determine, in each grid cell and 
at successive time intervals, thermodynamic characteristics of each non-aqueous phase 
and data representative of phase displacements in the reservoir ; and 
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- carrying out a delumping operation in order to obtain the detailed composition of the 
fluids in each cell. 

In order to characterize the thermodynamic path, the variations in the composition 
of the fluids and in the properties of the various phases are for example determined. 

5 During B.O. type simulation, thermodynamic parameters are for example 

determined in each cell, such as pressure (p™) , temperature T™) if it varies, the 
saturations of the liquid (SoJ) and vapour Sg") hydrocarbon phases, the injection or 
production rates, and, for each pair of cells (j,h), the volume flow rates of the liquid 
m™ a ) and vapour u™ h ) phases. 

10 The delumping operation comprises for example determining the equilibrium 

constants from input data specific to the delumping operation, and converting results 
expressed in volume into results that can be used in molar or mass quantities 
conservation equations. 

BRIEF DESCRIPTION OF THE FIGURES 

15 Other features and advantages of the method according to the invention will be 

clear from reading the description hereafter, given by way of non limitative example, 
with reference to the accompanying drawings wherein : 

- Figures 1 to 3 respectively show the variations, as a function of pressure, of three 
coefficients AD 0 , AD,, AD 2 during the differential vaporization of a reservoir oil 
20 simulated by means of Peng-Robinson's two-parameter equation. 
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DETAILED DESCRIPTION 

The various stages of the modelling method allowing notably to predict, as a 
function of time, the detailed composition of fluids produced by a hydrocarbon reservoir 
under production are as follows : 

5 a) The reservoir is first represented in form of a network of grid cells (j), each cell 

forming an elementary volume filled with fluid(s) in form of one or more phases, with 
at least one non-aqueous phase. The non-aqueous phases are referred to as hydrocarbon 
phases although they can contain certain components other than hydrocarbons, such as 
nitrogen, carbon dioxide, sulfur dioxide. 

10 b) The thermodynamic behaviour of the non-aqueous phase(s) required for a B.O. 

type thermodynamic representation is established in a way known in the art. The 
purpose is to describe the dependence of phase parameters on the pressure, the 
temperature if it varies during simulation, and possibly, for non-strict B.O. type 
modelling, an indicator of the evolution of the composition (for example the content of 

15 condensate dissolved in the gas phase) during a thermodynamic path (evolution) 
representative of the path that will be followed by the fluid in the reservoir. The 
commonest thermodynamic path is that of the differential operation during which the 
pressure of the mixture is progressively lowered at constant temperature : below the 
saturation point, one of the phases is removed at each pressure stage (for the differential 

20 vaporization of a liquid phase, the gas phase is removed) ; during this operation, the 
compositions and the properties of the phases (viscosities, densities, volume factor, gas- 
oil dissolution ratio in the liquid phase, etc.) are detennined. 
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It can be noted that it is also possible to deal with cases that require the use of 
several Black Oil thermodynamic representations for example, if the composition at the 
beginning of the simulation varies with the depth or laterally within the reservoir, or if 
several local thermodynamic paths can be distinguished during modelling. Several 
5 zones or ranges of variation of the thermodynamic or compositional quantities, often 
referred to as thermodynamic zones by specialists, can thus be defined and used. 

The origin of the data relative to the behaviour of the reservoir fluid can be of 
experimental nature (laboratory experiments), it can result from measurements in the 
field, or predicted from raw data and from complementary hypotheses by numerically 
10 simulating the behaviour of the reservoir fluid. Examples of the phase properties thus 
usually described as input data for B.O. type modelling are the viscosity of the phases, 
the gas-oil dissolution ratio, the volume factor of the oil, the compressibility factor of 
the gas. 

c) The fluid is defined, prior to B. O. simulation and for each thermodynamic zone 
1 5 or range, by a detailed representation with N«i components and/or pseudo-components. 

d) For each thermodynamic zone, an equation of state allowing to first reproduce, 
prior to B. O. simulation, the thermodynamic behaviour of the fluid during the 
thermodynamic path representative of the path that will be followed by the fluid in the 
reservoir is determined. 

20 e) The stages (often pressure stages) of the thermodynamic path of each zone are 

reproduced with the equation of state, in particular below the saturation pressure where 
two phases are present. This prior thermodynamic simulation uses, at the saturation 
pressure of the fluid and below, at each stage of the thermodynamic path, the 
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equilibrium constants obtained from the simulation results in order to calculate 
coefficient AD 0 and the n coefficients AD P of Equation (1) corresponding to the equation 
of state selected to describe the behaviour of the hydrocarbon fluid for example by 
mmirnizing the function : 


(2) O(AD 0 ,AD x ,...,AD p ,...,AD n ) = £ 


5 Equation (1) being obtained from an equation on the fugacity of constituent i in a 

phase where coefficients D 0 and the n coefficients D p specific to the phase appear : 

It is possible to alternately determine these coefficients that appear as phase 
parameters, the drawback being that the volume of input data is increased. 

We thus have, at each stage of the thermodynamic path considered to be 
10 representative of the path that will be followed by the fluid in the reservoir, the 
dependence as a function of the pressure, of the temperature, and possibly other 
necessary indicators, the properties of the phases necessary for B.O. simulation and, 
below the saturation pressure, coefficient AD 0 and the n coefficients AD P (or the n+1 
parameters D„ and D p , p ranging from 1 to n, by phase) that will be used to calculate the 
15 equilibrium constants during the delumping stage described hereafter. 

The variations of differences AD P (p = 0 to 2) for a differential operation on a 
reservoir oil simulated by means of Peng-Robinson's two-parameter equation are 
illustrated in Figures 1 to 3. It can be seen that the behaviour of these parameters is 
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sufficiently continuous to be readily introduced in form of a chart, or by a simple 
correlation as a function of the pressure. For example, in the present case for AD 0 : 

AD 0 = 1,406231 (1 -LN(P/Psat) [1 + 0,038 x LN(P/Psat)]). 

It can be convenient to add to the data stored for later use relative to the evolution 
5 of the values of these coefficients as a function of the pressure, temperature and 
compositional indicator abscissas, the evolution of the values of the oil density in the 
thermodynamic conditions of each stage of the thermodynamic path followed, above as 
well as below the saturation pressure. It is also possible to store, in addition to or instead 
of this data, other parameters of the oil and/or gas phases, but this is not essential. 

10 f) The thermodynamic behaviour of the hydrocarbon phase(s) is transposed, in a 

way known in the art, into one or more B.O. type thermodynamic representations) in 
the input data of the model. These dependences are described above the saturation 
pressure of the mixture of hydrocarbons and below the saturation pressure, either point 
by point in charts with an interpolation method and, if need be, an extrapolation method, 

15 or by analytical correlations, or by a combination of charts and correlations. 

g) A B.O. type simulation is carried out in a way known in the art, said simulation 
allowing to compute at least in each grid cell (j) and at successive time intervals (m, 
m+1, etc.), pressure (p™) , temperature (T™) (if it varies), the liquid (SO™) and vapour 
(Sgj) hydrocarbon phase saturations, the injection or production rates, and, for each 

20 pair of cells (j,h), the volume flow rates of the liquid (u™ h ) and vapour (u" h ) phases. 

h) The equilibrium constants (k™) of component i are determined, at each time 
interval (m) and in each cell (j), from the dependences of coefficients (AD 0 , AD b ..., 
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AD P ,..., AD n ) of Equation (1), or of coefficients (D 0 , D 1; ..., D p ,..., D n ) of Equation (3) 

by phase, established at e) with pressure (/>*) , temperature (Tf) and, if need be with 

non-strict B.O. modelling, with a composition indicator such as the content of 
condensate dissolved in the gas phase. 

i) The vaporized fraction (<?") is determined at each time interval (m), in each cell 
(j), either from the B.O. type simulation results, or, for higher accuracy in case of 
coexistence of the oil and gas phases, by solving Rachford-Rice's equation known to 
specialists, which is applied for example in the following document : 

- Rachford H.H.Jr and Rice J.D. ; « Procedure for Use of Electronic Digital 
Computers in Calculating Flash Vaporization Hydrocarbon Equilibrium », J. Pet. 
Technol., 1952, 14, 19, 

from the molar fractions of each component i in the global detailed composition (z™) of 
the hydrocarbon fluid in cell (j) at time interval (m) 


j) The detailed composition of each hydrocarbon phase is estimated at each time 
interval (m) and in each cell (j) by means of the following relations : 


(4) Z 


= 0. 


(5) *,; = 


for the oil phase 


i+(*;-i)0; 


(6) y ; = 


for the gas phase. 
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k) The molar density of the oil phase E,o m } is estimated in cell (j) at time interval 
(m) for example from the density of the oil po™ , in cell (j) at time interval (m), from 
the following relation : 


" J MMo 


(7) $>; = - 

MMo being the molar mass of the liquid hydrocarbon phase (referred to as « oil »). 

In Equation (7), density po™ can therefore be obtained from the results, at time 
interval (m), of the B.O. simulation or from the input data specific to the delumping 
stage and stored during the prior thermodynamic simulation. MMo can be calculated 
from the results of the B.O. simulation if the molar mass of the heavy hydrocarbon 
constituent is defined by default or in the user input data, or from : 

Nrd 

(8) mm?; = £mm ; ^ 

MM; being the molar mass of component or pseudo-component i, and quantities x™ 
being calculated in the previous stage. 

1) The molar density of the gas phase is estimated in each cell (j) at each time 
interval (m), either from input data specific to the delumping stage and stored during the 
prior thermodynamic simulation, or from the results of the B.O. simulation at pressure 
PP , temperature r™ , for example, when the compressibility factor of the gas is 

included in the input data of the B.O. type thermodynamic representation, with the 
following relation : 
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< 9) *>*gif 

Zg™ being the compressibility factor of the gas, R the perfect gas constant, 

or, for example, if it is the volume factor of the gas that is included in the input data of 
the B.O. type thermodynamic representation, with the following relation : 

(io) &; = - 1 


jVmolst 

Vmolst being the molar volume of the gas under standard conditions and Bg™ being the 
5 volume factor of the gas. 

m) The molar flow rates of the liquid (u™ h ) and vapour (w" A ) phases are evaluated 
in each cell (j) at each time interval (m) from the volume flow rates of the liquid (u" k ) 
and vapour (u™ h ) phases obtained during B.O. simulation and the molar densities 
obtained in stages k) and 1). 
10 n) The molar fraction of each component or pseudo -component i is estimated at 

time interval m+1 in the global detailed composition (z ; J +1 ) of the hydrocarbon fluid in 
cell 0"), knowing the molar fraction (z") at time interval m, from the following 
equations : 

z;n; - utys; + x ;s;) - At £ + ) 

(1 1) z m+1 = ^Jl 

(12) N? 1 = N™ - At(S™ +S™)-At 5>" +«" ) 
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Equation (12) expressing the total hydrocarbon molar balance on cell (j), iV™ being the 
hydrocarbon molar amount contained in cell (j) at time interval (m), taking into account 
the matter exchanges with all the cells (h) adjoining (j) which form set J(j). 

In Equation (11), the writing of terms (y™) and (jc™) - wherein j'=j for a flow from 

5 cell (j) to cell (h) or in the well, and j'=h for a flow from cell (h) to cell (j), j' 
corresponding to the fluid injected in the case of injection wells, S being then negative - 
implicitly presupposes the use of a simple upstream pattern for the compositional flows. 
A more general writing of these terms is (y™ h ) and (x™ h ), where x™ h and y™ h describe 
the compositions of the liquid and gas phases, obtained in a way known in the art, 
10 flowing between cells (j) and (h). 

Source terms (molar flow rates) (S™) and (S™) are obtained from source terms 
(volume flow rates) (S™) and (5") in the B.O. simulation by repeating stages j) to m) 
for the source terms. 

Stages h) to n) describe the delumping operation. The details provided are based on 
15 equations of conservation in moles ; they could just as well be based on equations of 
conservation in mass. 

From the moment that one knows how to describe the detailed composition in each 
cell of the B.O. type simulation model at the time t, which can be, in particular, the 
initial time, one can describe, by means of the present method, the evolution of the 
20 detailed composition in each cell during the development process modelled in a « Black 
Oil » type model. 


